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P3 \ ing different standard phonon basis obtained through the Lang Firsov (LF), modified 

O ' 

^ ' LF (MLF) and modified LF transformation with squeezed phonon states (MLFS). 

^ , The role of these convergent expansions using the above prescriptions in lowering 

\ the energy and in determining the correlation functions are compared for different 

values of coupling strength. The single-electron energy, oscillator wave functions and 
correlation functions are calculated for the same system. The applicability of differ- 
ent phonon basis in different regimes of the coupling strength as well as in different 
regimes of hopping are also discussed. 
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The two-site single-polaron problem is studied within the perturbative expansions us- 
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1. Introduction 



The Holstein model [0 is one of tiie fundamental models describing the inter- 
actions of conduction electrons with lattice vibrations. The model in the simplest 
form consists of one electron hopping term, dispersionless phonons and an interaction 
term which couples the electron density and ionic displacements at a given site. The 
interaction term favours localization of the electron. When the interaction is strong, 
the gain in the localization energy outweighs the kinetic energy of the electron and a 
self-trapping of the electron (or a hole as the case may be) occurs with the creation of 
lattice deformation in the immediate vicinity of the charge carrier. The motion of the 
electron is then accompanied by the lattice deformation which results in a reduced 
effective hopping of the dressed electron (polaron). In the antiadiabatic limit, where 
the phonon frequency (cuq) is greater than the electronic hopping integral (t), the con- 
finement of the lattice deformation around the charge carrier is very local for large 
electron-phonon (e-ph) interaction. This gives rise to small polarons whose nature 
and dynamics is generally studied using the Lang-Firsov (LF) method based on the 
canonical LF transformation 0. Recently, Ranninger and Thibblin ^ and de Mello 
and Ranninger studied a two-site Holstein model using the numerical diagonal- 
ization technique and obtained non trivial results which are difficult to understand 
from the classical zero phonon averaging LF approach. In fact, in Ref. the au- 
thors expressed doubts regarding the applicability of the LF approach and validity of 
the strong coupling perturbation expansion (in hopping) even in the strong coupling 
antiadiabatic limit where it is believed to be valid. Considering a similar two-site sys- 
tem Firsov and Kudinov studied the energy levels and the wave functions of the 
system within the perturbation theory and concluded that the results are consistent 
with the exact numerical results obtained in Refs. and Q for large e-ph coupling 
strength and small hopping. Marsiglio ^ studied the Holstein model in one dimen- 
sion with one electron up to 16 site lattices using numerical diagonalization technique 
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and concluded that neither the Migdal nor the usual small-polaron approximation 
is in quantitative agreement with the exact results for intermediate coupling strength. 
Kabanov and Ray |^ and Alexandrov et al. performed the exact numerical diag- 
onalization of two, four and six sites with one electron and observed that for t > ujq 
the adiabatic Holstein small-polaron approximation, whereas for t < ujq the LF ap- 
proach describes the ground state energy of the system accurately except for a region 
of intermediate coupling strength. So no single conventional analytical method at 
present is known to us so as to describe the Holstein model for the entire range of 
the coupling strength either in the adiabatic or in the antiadiabatic limit. In a recent 
work we addressed this problem and considering a two-site system we developed 
a perturbation expansion using modified LF phonon basis. We calculated the per- 
turbation corrections to the ground state energy and wave functio up to the third 
and second order respectively. The results obtained thereby are in good agreement 
with the exact numerical results for the entire range of the coupling strength for the 
intermediate value of hopping (t ~ t^o). 

In this paper we have extended our previous work by calculating up to fifth order 
perturbation correction to the wave function and sixth order correction to the energy 
for the ground state of the two-site single polaron system using various standard 
phonon basis obtained through the LF, modified LF (MLF) and modified LF with 
squeezing transformations (MLFS). It may be mentioned that the role of squeezing is 
believed to be important for intermediate coupling and hopping regime |TT]-|T^ and 
studies on a two-site and a four-site Holstein model using the MLFS transformation 



followed by zero phonon averaging show that the energy obtained within such 

method is very close to the exact result. 

In the LF approach a phonon basis of fixed displacement at the electron residing 
site is chosen. Such a choice of basis diagonalizes the Holstein hamiltonian in absence 
of hopping. The hopping term is then treated as a perturbation . The perturbation 
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series is naturally expected to converge when hopping is weak and the e-ph coupling 
is strong. However, it is not precisely known the limit of the coupling strength as a 
function of hopping beyond which the LF approach is valid. Within the LF approach 
in the zeroth order of perturbation the effect of retardation between the electron 
and the lattice deformation produced by the electron cannot be obtained. Such a 
retardation becomes very important with smaller e-ph coupling strength and larger 
hopping. The phonon basis chosen in the MLF or MLFS approach have variational 
displacements and can produce the retardation effect even in the zeroth order of 
perturbation [|T^,n|. The ground state energy predicted within the MLF and MLFS 
method in the zeroth order of perturbation is much lower than that within the LF 
approach except for large values of the coupling strength, where all the methods 
become equivalent. All these suggest that the MLF or MLFS phonon wave functions 
are proper choices for perturbative calculation when hopping is appreciable and the 
coupling strength is not very high. For weak hopping and large values of the coupling 
strength the MLF or MLFS phonon wave function reduces to that of the LF approach. 
Thus the perturbation methods based on the MLF and the MLFS phonon wave 
functions have the potentiality to be applicable for a wide range of the coupling 
strength. 

The objective of this work is to develop perturbative expansions for the two-site 
single polaron system using phonon basis obtained through the LF, MLF and MLFS 
transformations and investigate the convergence of the perturbation series for different 
physical quantities of interest as a function of the e-ph coupling strength and hopping. 
Such study would provide direct answers to the important questions as (i) what is 
the region of e-ph coupling strength where the LF method is valid, (ii) whether the 
MLF and MLFS phonon wave functions are better than the LF wave function in 
the zeroth order of perturbation, (iii) in which region MLF and MLFS methods give 
better results than the LF method, (iv) which method within a few orders (one or 
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two) of perturbation could predict results closer to the exact one. 

In Sec. II we define the model hamiltonian and describe different variational 
phonon basis states that we have considered. In sec. Ill we present the expressions 
for the energy, wave function and static correlation functions calculated within the 
perturbation theory for the ground state. In Sec. IV we shall present the results 
obtained by different methods and discuss about the convergence of the perturbation 
series, hence the applicability of the methods in different regions of the e-ph coupling 
strength for different values of hopping. 

II. Formalism 

The two-site single-polaron Hamiltonian is 

where i =1 or 2, denotes the site. Qo- (cJ^.) is the annihilation (creation) operator 
for the electron with spin a at site i and riia (=c|^Cio-) is the corresponding number 
operator, g denotes the on-site e-ph coupling strength, t is the usual hopping integral. 
hi and h\ are the annihilation and creation operators, respectively, for the phonons 
corresponding to interatomic vibrations at site i, ujq is the phonon frequency. The 
hamiltonian (1) has spin degeneracy, hence for the one electron case the spin index 
is redundant. So spin indices are dropped in the following. 

Introducing new phonon operators a — (6i -|- h-})l\pl and d — {bi — &2)/'\/2 the 
Hamiltonian is separated into two parts {H — Hd-\- Hg) : 

Hd = Y^ erii - t{c{c2 + 4ci) + uJog+{ni - n2){d + d^) + UodU (2) 

i 

and 

Ha = ouoo'd - ujQ-n?g\ (3) 
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where g+ = g/ -\/2 , 5 = a + ng+ and a) = + ng+ . 

Ha describes a shifted oscillator which couples only with the total number of 
electrons n{= rii + 77-2), which is a constant of motion. The last term in Eq.(3) 
represents lowering of energy achieved through the lattice deformations of sites 1 and 
2 by the total number of electrons. 

Hd represents an effective e-ph system where phonons directly couple with the 
electronic degrees of freedom and its solution by any analytical method is a non 
trivial problem. We now use the MLF transformation where the lattice deformations 



produced by the electron are treated as variational parameters ||r^ , P~4|JT7[| . For the 
present system, 

Hd = e^i/^e-^ (4) 

where R = X{ni—n2){S — d), A is a variational parameter related to the displacement 
of the d oscillator. 

The transformed Hamiltonian is then obtained as 

Hd = ujod'^d + ^ epTii - t[c[c2 exp(2A(d'^ - d)) 

i 

+ 4ci exp(-2A(dt - d))] + ^o(g+ - A)(ni - n2)(d + d^) (5) 

where 

ep = e - cuo(2fl'+ - A)A (6) 

For a perturbation method it is desirable to use a basis where the major part of 
the hamiltonian becomes diagonal. When the hopping is appreciable a retardation 
between the electron and associated lattice distortion sets in and as mentioned before 
the MLF or the MLFS method would work better than the LF method for a wide 
region of e-ph coupling strength. 

Now we will make a squeezing transformation |18] to the Hamiltonian (5) 



Hsd = e'^HdC-^ (7) 



where S = a{didi — d\d\). The new phonon basis is squeezed with respect to the 
previous basis. Squeezing parameter (a) partly reduces the polaronic narrowing effect 
and consequently enhances the hopping. However, the phonon energy increases with 
increasing a and a competition between phonon energy and hopping delocalization 
(kinetic) energy determines the value of a The transformed hamiltonian (7) 

takes the form 

Hsd = uJodU[cosh^{2a) + sinh^(2a)] + ujq cosh(2a) sinh(2a)(dd + dM'^) 
+ E ^p^* - t[c\c2 exp(2Ac(dt - d)) + c^Ci exp(-2Ae(dt - d))] 

i 

+ ujQ{g^ — \){ni — n2){d + d^) exp(2Q;) + ujq sinh^(2a) (8) 

where Ac = A e"^" (9) 
For the single polaron problem we choose the basis set 

l±,A^) = ^(cI±4)|0)e|iV) (10) 

where |+) and |— ) are the bonding and antibonding electronic states and |A^) denotes 
the iVth excited oscillator state in the MLFS, MLF or LF basis depending on the 
method considered. It may be noted that the MLFS basis turns into the MLF basis 
if one puts a = and the MLF into the LF basis when A = (7+. 

The hopping term Ht in Eq. (8) has both diagonal and nondiagonal matrix 
elements in the chosen basis (10). The diagonal part of Ht is given by 

{N,±\Ht\±,N) = {N,±\-t[c\c2 exp(2Ae(d"^ - d)) + c^ci exp(-2Ae(d"^ - d))]|±, N) 



N 



-1)W< 



i=0 

„„j Ar_ _ N\ 
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where te = t exp(-2AJ and Na, - i,(;v-i)!- 

For the perturbation method we consider the diagonal part of the Hamiltonian 
(8) as the unperturbed Hamiltonian (Hq) and the remaining part of the Hamiltonian, 
Hi = Hgd — Hq, as a perturbation. 



The unperturbed energy of the state |±, -/V) is given by 



£5v = ±l^o|±, N) = a;o[sinh2(2a) + N(sinh2(2a) + cosh2(2a))] 

.i=0 



(12) 



The general off-diagonal matrix elements of Hi between the two states |±,-/V) and 
|±, M) are calculated as (for {N - M) > 0) 



(TV, ±\Hi\±,M)= P{N, M) + y ^iV(7V - 1) sinh(4a) 5n,m+2 (13) 

for even (N - M). 



(TV, ±|//i|t, M) = P{N, M) + VNuj^{g+ - \)e'"'5N,M+r 

for odd (N-M). 



where P(N, M) = ^ie(2Ae) 



N-M , 



M! 



M 



{N-M)\ 



+ E[(-i)'' 



R=l 



^^'^"''''^ -M(M-1)...(M-R + 1)] 



{N - M + ny.Ri' 



(14) 



III. Perturbation corrections to the ground state 



In this section we present the calculations of the perturbation corrections to the 
ground-state energy and wave function up to the sixth and fifth order, respectively 
using the LF, MLF and MLFS methods. The static correlation functions relating the 
electron and associated lattice deformations are also calculated using the correspond- 
ing perturbed wavefuctions. 

For the chosen basis (10) the state |+)|0) has the lowest unperturbed energy, 

Eq'^ = ep — te + loq sinh^(2Q;). The matrix element connecting this ground state and 

an excited state |e, N) is given by 

(2A,^^ 



(iV,e|ifi|+,0) 



+ LJo{9+ - A)e " 6n,i 



for odd N 



(15) 



(2X ) 

-te—y=- + V^Uo sinh(2Q;) cosh(2Q;) 5n,2 



5e,+ for even N 
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The first order correction to the ground state wave function is obtained as, 
I ,(1), _ [cuo((7+ - A)e2" - 2Aete] I , [-te(2Ae)' + 2cuo sinh(2a) cosh(2a)] 



^e(2Ae'*^ 



where -E'^^^r is the unperturbed energy of the state |±, A^) as given in Eq.(12) and e 
= + or - for even and odd N respectively. 

The first order correction to the energy {Eq'^) is zero since Hi has no diagonal 
matrix element. The second order correction to the ground state energy is given by 

—jm ^ ^o{9+ - Aj e dN,l\ 



(2) _ I I 



Af=l,3.. l-C/Q - -C'-.Arj 

I^M^ + 72 cuo sinh(2a) cosh(2a) 5n,2|' 
+ E ,^(0) ^(0). (17) 

iV=2,4.. l-C/Q - -^+,ArJ 

Now one has to make a proper choice of A and a so that the perturbative expan- 
sion becomes convergent. Usually, within the MLF and MLFS approach zero phonon 



averaging is made |12,14] and the variational parameters A and a are found out by 
minimizing the ground state energy of the system. This corresponds to the minimiza- 
tion of the unperturbed energy in our calculation. Previously we have followed this 



procedure within the MLF approach \U]\ and found that the perturbation corrections 
to the energy converge rapidly. Here we will follow the same procedure to find out 
the variational phonon basis as a function of e-ph coupling for the MLF and MLFS 
methods. Minimizing the unperturbed ground state energy Eq'^ with respect to A 
and a we obtain 



tuo + 2tee-4" 
a = isinh-^[2a(g+-A)] 

It is interesting to note that these choices of A and a make the offdiagonal matrix 
elements between the states |+, 0) to |— , 1) and |+, 2) equal to zero and consequently. 



the coefficients of |— , 1) and |+, 2) in Eq.(16) and tlie terms corresponding to = 1, 2 
in tfie r.fi.s of Eq.(17) vanisfi. Since tfie first two terms in Eq. (16) are zero for 
tlie MLFS choice of the basis it is expected tliat the first order correction to the 
ground state wave function and second order correction (Eq. 17) to the energy would 
be smaller within the MLFS method. To examine the validity of the perturbation 
procedure one should calculate the higher order corrections and check whether the 
series is converging properly. 

The second order corrections to the ground state wave function is given by 



{2)\ 



E4?|k) 

(Hi)km(Hi)mO T 



where ag,-* = ^ [ 



(18) 



(0)^^17(0) 



where |k)'s denote the unperturbed states |±, A^) with the unperturbed energy E^^^ 



±,N 



and the refers to the ground state |+, 0). (i?i)km is the off-diagonal matrix element 
of Hi between the states |k) and |m). These matrix elements are given in Eqs. (13) 
and (14). The third order correction (Eq^^) to the ground state energy is given by, 



4'^ = E(^i)ok4? 



(19) 



In general the nth order correction to the wave function and the (n-l-l)th order 
correction to the energy are given by 



fc^O 

^(n+i) ^ ^(^Hi)okaQl^ 



(20) 
(21) 



where 



E, 



(0)^ 



_p(l)„(n-l) 
-^0 %]si 



-^0 ^Ok 



p(3) {n-3) 
-^0 %k 



(n-l)n 



-^0 ^Ok 



(22) 
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Using Eqs. (20) and (21) all the higher order corrections to the wave function and 
energy may be calculated step by step. 

The ground state wave function may be written as, 

|V'G)^l+,o) + |#) + |V'f ) + !#) + 

= |+,0)+ Yl cn\-,N)+ Y1 (23) 

iV=l,3,.. N=2,4,.. 



The coefficients cjv are determined from the sum of the corresponding coefficients a'^^, 
«ok > «ok> etc., wh( 
function \iPg)n is 



^ok ) ^ok ) etc., where N represents the state |k). The normahzed ground state wave 



|V'G)iV = 



1 



where Nq is obtained as 



G 



7V=1,2,3,.. 



Correlation function 



The static correlation functions {niUi)o and (711^2)0, where Ui and U2 are the 
lattice deformations at sites 1 and 2 respectively, produced by an electron at site 1, 
are the standard measure of polaronic character and indicate the strength of polaron 
induced lattice deformations and their spread. The operators involving the correlation 
functions may be written as 



niu^,2 = + ±id + d'^y - 2 {ng+ ± A(ni - 712))] 



The final form of the correlation functions are obtained as 

(niMi)o = ^ -{g+ + A) + 

(niM2)o = 7; -{9+ - A) - 



where Aq = (^/^G|ni(d + d'^)|t/^G) = E "^n \/N cn-i + VN + 1 cn+i 



A, 


g2a- 


Ng 


A, 


g2a- 


Ng _ 


Cn 





(24) 



(25) 
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IV. Results and discussions 



In this paper we present the results for t—0.5, 1.1 and 2.1 (in a scale of cuq =1). 
These values of {t/uio) covers the cross-over region from the antiadiabatic to the 
adiabatic limits. Furthermore, exact results on some of the physical quantities are 
available for the above values of {t/uo) which enables us to compare our results 
obtained within the perturbation methods with the exact results. 

We estimate the perturbation corrections to the ground state energy and wave 
function up to the sixth and fifth orders respectively, within the LF, MLF and MLFS 
methods considering 25 phonon states (which is sufficient for < 2.2) in the trans- 
formed phonon basis. Calculations of such successive orders of corrections would 
provide us the key features of the convergence of the perturbation series and the 
applicability of the method concerned. 

For (t/cuo) — 0.5 we find a rapid convergence in the perturbation series for the 
ground state energy for all values of g+. Fig. 1 shows the relative perturbation cor- 
rections i.e., the ratios of the perturbation corrections (to the ground state energy) of 
different orders to the unperturbed energy as a function of g+. The energy corrections 
within the MLFS and the MLF methods are much smaller than that within the LF 
method for < 1.6. For {t/uJo) = 0.5 the fifth or sixth order correction to the energy 
is so small that the energy obtained with fifth or sixth order of perturbation within 
MLF method may be treated as the exact energy for all values of g^. In Fig. 2(a) we 
plot the exact energy thus obtained along with the unperturbed energies within dif- 
ferent methods. It is seen that the MLFS or the MLF unperturbed energies are closer 
to the exact one than the LF unperturbed energy. In Fig. 2(b) we show the variation 
of |('?/'o°'*I'V'g)P) which gives the probability that the perturbed ground state {\iIjg)) lies 
in the unperturbed ground state \^pQ^^), with g+ (IV'g) is found out considering up to 
the fifth order correction to the wave function) . This probability is closer to 1 for the 
MLFS and the MLF methods compared to that in the LF method. The results of Fig. 
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2 clearly indicate that the phonon basis, chosen within the MLFS or the MLF meth- 
ods, are much better than the LF basis. It may be mentioned that the correlation 
function {niU2)o, which is a measure of the lattice deformation at site 2 produced by 
an electron at site 1 due to retardation effect vanishes for the LF unperturbed state. 
This correlation function is very sensitive to the higher order perturbation corrections 
to the wave function than the ground state energy. To examine the convergence of 
the perturbative expansion on this correlation function we plot {niU2)o, calculated 
up to different orders of perturbation corrections, as a function of g+ in Fig. 3. It 
is found that for (t/o^o) = 0.5 this correlation function converges for all values of g+. 
The curves corresponding to the correlation function calculated up to fourth and fifth 
order corrections to the wave function almost merge together (Fig. 3) and specially 
for the MLF method they are almost identical in the entire region of g^. Hence, the 
correlation function obtained considering up to the fifth order correction within the 
MLF method, may be considered as the exact one for {t/coo) — 0.5. Fig. 3 also shows 
that even in the strong coupling limit one has to consider up to the second order 
corrections to the wave function to obtain (nearly) exact results. It is interesting to 
note that the MLF and the MLFS methods can predict non zero values of (711X^2)0 
even in the zeroth order of perturbation (i.e. with the unperturbed ground state wave 
function in the MLF or MLFS basis) for low and intermediate values of g+ whereas 
the unperturbed LF wave function always predicts zero value for (7111*2)0 and we find 
that in an appreciable region of low values of g^ the predictions of the MLFS method 
with zeroth order of perturbation is very close to the exact one. 

In Figs. 4-8 we have shown the results for t/ujQ — 1.1 which is in between the 
adiabatic and antiadiabatic limits. In Fig. 4 the relative perturbation corrections to 
the unperturbed energy are plotted as a function of g+. The perturbation corrections 
in energy are smaller within the MLF and MLFS than the LF method. For each 
method there is a range of g+ where higher order corrections are appreciable and the 
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convergence is weak. For t/uoQ = 1.1 the range of (7+ where the convergence is weaker, 
is 0.4-1.2 for the LF, 1-1.3 for the MLF and 1.1-1.55 for the MLFS method. For 
values of g+ outside this range the perturbation corrections show rapid convergence. 
In the MLF method the perturbation corrections show satisfactory convergence in the 
entire region of (7+, while for the LF and MLFS methods there are small regions of 
where fifth or sixth order corrections are comparable or greater than the third and 
fourth order corrections and the convergence is not satisfactory in that region. The 
ground state energy calculated within the MLF method considering up to the sixth 
order correction, is shown in Fig. 5 as a function of g+. This energy when compared 
with the exact results of Ref. is found to be almost identical with the exact 
one. This is expected because the higher order corrections are successively smaller 
in magnitude and of alternate sign within the MLF method, so the net contributions 
from seventh to higher order terms to the energy would be negligibly small with 
respect to the unperturbed energy. The unperturbed ground state energy for the 
LF, MLF and the MLFS methods are also shown in the same figure. The MLFS 
unperturbed energy is closer to the exact energy and for low values of it is almost 
identical to the exact energy. 

In Fig. 6 we plot the correlation function {niU2)o obtained by considering up 
to the different orders of perturbation corrections to the wave function against (7+. 
The MLFS method shows excellent convergence for low values of (7+ (< 0.9), the LF 
method shows very good convergence beyond (7+ = 1.2. The MLF method shows good 
convergence for the entire region of g^ , however the convergence is weaker in the range 
0-9 < g+ < 1.3. When compared with exact results of {niU2)o (taken from the Ref. 
1^), it is found that the MLF results up to the fifth order perturbation are identical 
to the exact results except in the region 0.9 < g^ < 1.3 where a slight departure in 
values from the exact results is seen. The correlation function {niUi)o would evidently 
show a much more rapider convergence than the convergence of (niM2)o- 
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To examine the nature of the cross-over from the delocahzed (large) to locahzed 
(small) polaron we plot {ni{ui — U2))o/g+ against g+ in Fig. 7 where the (nearly) 
exact plot (obtained from the MLF considering up to fifth order corrections to the 
wave function) and those obtained from the MLF and MLFS in the zeroth order of 
perturbations are shown. The standard LF approximation (zeroth order of perturba- 
tion or zero phonon averaging) always predicts the value of this quantity to be equal 
to 1, which is a characteristic of extremely localized polarons, and so a cross-over 
from small to large polaron cannot be obtained in the standard LF approximation. 
The MLF and MLFS methods can predict this cross-over even in the zeroth order 
of perturbation [|r^,0,|r^. It may be mentioned that the value of {ni{ui — U2))o/g+ 
within MLF and MLFS methods in the zeroth order of perturbation is X/g+ (Eq. 
(25)). Fig. 7 shows that up to g^ = 1.0 the results of MLFS with zero phonon 
averaging almost coincides with the exact results. The (nearly) exact plot shows 
a smooth cross-over from the delocalized to localized polaron with increase of g^, 
whereas the zero phonon averaging results within MLF or MLFS procedure shows an 
abrupt feature which were noted before |12||1^,|15| . Lowen pointed out that for a 
finite phonon frequency there cannot be any abrupt cross-over in the ground state of 
an e-ph system. The (nearly) exact plot of Fig. 7 is consistent with this conclusion. 

In Fig. 8 we have shown the ground state wave function (within the MLF method 
considering up to fifth order correction to the wave function) for the d oscillator as 
a function of position x for different values of the e-ph coupling when the electron is 
located at site 1. For weak coupling the wave function shows displaced Gaussian like 
single peak. For intermediate values of g^ an additional prominent shoulder appears 
in the wave function as seen in the curve for g^=1.3. For higher values of (?+ this 
shoulder takes the form of a broad peak. These results are completely consistent with 
the results obtained by Ranninger and Thibblin by exact diagonalization study 0. 

In Figs. 9-12 we have given the results of t/cuo = 2.1, which is towards the adiabatic 
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limit where the LF method is not expected to work in this region. So it would be 
of general interest to examine whether the methods, dealt in this paper, show good 
convergence in any region of In Fig. 9 the relative perturbation corrections of 
different orders to the ground state energy are shown as a function of g+. The figure 
shows that the perturbation corrections within the LF method are very large and 
do not converge at all in an appreciable region of (0.6-1.25). For the MLF and 
MLFS methods the perturbation corrections in energy are much smaller. The MLF 
method shows good convergence except in a region of 1.2 < g+ < 1.6. For the MLFS 
method the corrections are negligible in a wide region of 5^+ {0 < g+ < 1.1). However, 
in a region 1.25 < g^ < 1.7 the perturbation corrections within the MLFS do not 
converge. For g^ > 1.4 the convergence within the LF method is satisfactory and 
hence may be used. 

In Fig. 10 we show the variation of KV'o^^IV'g)!^ with g+ for different methods. 
The figure shows that for t/ujo — 2.1 this value within the LF method is far from 
1.0 for a broad region of g+ (0.4—1.3) which signifies that the LF method is totally 
inapplicable in that region. However, in the same region the value of \{ipQ^\ipG)\'^ 
within the MLFS is close to 1 and thus MLFS should work exceedingly well there. 

In Fig. 11 the convergence of the correlation function {niU2)o is shown. It is 
clearly found that the MLFS method shows excellent convergence for low values of 
g+ (< 1.2), the LF method shows good convergence beyond g+ — 1.5. The MLF 
method shows a systematic good convergence for the entire region of g+, however the 
convergence is weaker in an intermediate region of g^ . 

In Fig. 12 the ground state phonon wave functions, evaluated by considering 
up to the fifth order perturbation correction within the MLF and LF methods, are 
shown for the d oscillator as a function of position x for = 1.6 and for different 
values of t/coo (0.5, 1.1, 2.1), when the electron is located at site 1. The LF and 
the MLF predictions are identical and hence cannot be distinguished in the plot. 
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These plots also exactly match with the phonon wave functions, calculated with 
exact diagonalization technique 0], for the corresponding set of parameters. The 
reason for such excellent agreement with the exact results is that for = 1.6 and 
for t/uo = 0.5, 1.1 and 2.1, a very good convergence is achieved in the wave function 
within the LF and MLF methods which is evident from Figs. (3), (6) and (11). 

V. Conclusions 

In the present work we examine the convergence of the perturbation expansions 
using the LF, MLF and MLFS phonon basis for a two-site one electron Holstein model 
for three values of t/uo (0.5, 1.1 and 2.1). For t/uQ = 0.5 the perturbation series 
show very good convergence within all the three methods in the entire region of e-ph 
coupling strength. The perturbation corrections within the MLF and MLFS methods 
are much smaller than that of the LF method. For higher values of t/ujQ also good 
convergence is achieved for the major region of g^. The perturbation expansion shows 
weak or bad convergence in an intermediate region of g+, the range of which varies 
in different methods. The MLFS method yields negligible perturbation corrections 
and shows good convergence for a wide region of (7+ (from zero to an intermediate 
value) which increases with increasing value of t/ujQ. In this region the LF method 
particularly shows bad convergence hence inapplicable. For large values of g^ the LF 
method shows good convergence, as expected; the MLF and the MLFS methods also 
work satisfactorily there. The MLF method shows a systematic convergence in the 
whole region of g+ even for higher values oft/uo, as considered here. It is interesting 
to note that for strong coupling strength consideration of the second order correction 
to the energy as well as to the wave function reproduces almost the exact results 
for all the three methods. The calculated oscillator wave functions following our 
method match excellently with the exact oscillator wave functions, derived from the 
numerical diagonalization studies. Our study also shows that the region of g^, where 
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no perturbation analytical method is applicable, is really very narrow in contrast to 
the conclusions of some earlier works I^Q. Although the analysis presented in this 
paper is based on a two-site system, the approach may permit a better progress for 
analytical studies of many-site systems. 
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Figure captions : 

FIG. 1. Variation of the relative perturbation corrections e'^^ /E^"^ to the ground 
state energy as a function of the couphng strength {g^) for t/uoQ = 0.5 for (a) LF, 
(b) MLF and (c) MLFS method. £^o"^ the nth order perturbation correction to the 
ground state energy and Eq*^ is the unperturbed ground state energy. 

FIG. 2. (a): The plot of the ground state energy obtained by considering pertur- 
bation correction up to fifth order within the MLF method as a function of for 
t/cuo — 0.5. The unperturbed energies (equivalent to the energy obtained with stan- 
dard zero phonon averaging (ZPA)) within the LF and MLF methods are also shown. 
The MLFS (ZPA) results (not shown) lies in between the MLF and the exact curves, 
(b): The variation of |(4°Vg)P with g+ within LF, MLF and MLFS methods. IV^g) 
and I ■00^^) a-re the corrected and the unperturbed wave functions, respectively. 

FIG. 3. Plot of the correlation function (^1^2)0 calculated up to different order of 
perturbations in the wave function versus g+ for t/utQ = 0.5 within different methods 

(a) LF, (b) MLF and (c) MLFS. The labels (2), (3), .... denote the curve obtained 
by considering up to the second, third, .... order correction to the wave function, 
respectively. 

FIG. 4. Variation of the relative perturbation corrections e'q'^ / Eq'^ to the ground 
state energy as a function of the couphng strength [g^) for t/ujQ — 1.1 for (a) LF, 

(b) MLF and (c) MLFS method, e'q'^ is the nth order perturbation correction to the 
ground state energy and £^0°'* the unperturbed ground state energy. 

FIG. 5. The plot of the ground state energy obtained by considering perturbation 
correction up to fifth order within the MLF method as a function ol g^ lor t / ujq — 1.1. 

The unperturbed energies (or with zero phonon averaging results (ZPA)) within LF, 
MLF and MLFS methods are also shown. 
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FIG. 6. Plot of the correlation function (72,1^2)0 calculated up to different order of 
perturbations in the wave function versus g+ ior t/cuo — 1.1 within different methods 

(a) LF, (b) MLF and (c) MLFS. The labels (2), (3), .... denote the curve obtained 
by considering up to the second, third,... order correction to the wave function, 
respectively. 

FIG. 7. The variations of {ni{ui - U2))o/g+ with g+ for LF, MLF and MLFS 
methods with zeroth order of perturbation for t/uio = l-l- The solid curve corresponds 
to that obtained within the MLF method up to fifth order correction to wave function. 

FIG. 8. Ground state oscillator wave function G{x) as a function of x for different 
values of the coupling strength when the electron is located at site 1 for t/cuo — 1.1. 

FIG. 9. Variation of the relative perturbation corrections E^"''' /Eq'^ to the ground 
state energy as a function of the coupling strength for t/ujo — 2.1 for (a) LF, 

(b) MLF and (c) MLFS method. E'q"^ is the nth order perturbation correction to the 
ground state energy and Eq'^ is the unperturbed ground state energy. 

FIG. 10. The variation of | (4°^ |^/'g) P with g+ within LF, MLF and MLFS methods 
ior t/uJo — 2.1. \i/jq) and |'0o°^) corrected and the unperturbed wave functions, 

respectively. 

FIG. 11. Plot of the correlation function (niU2)o calculated up to different order of 
perturbations in the wave function versus gj^ ioi t/ujQ = 2.1 within different methods 
(a) LF, (b) MLF and (c) MLFS. The labels (2), (3), .... denote the curve obtained 
by considering up to the second, third, .... order correction to the wave function, 
respectively. 

FIG. 12. Ground state oscillator wave function G{x) as a function of x for g^ = 1.6 
when the electron is located at site 1 with different values oit/ujQ. 



22 







LF (ZPA) 

MLF (ZPA) 
MLFS (ZPA) 




FIG. 6 




FIG. 9 



Authors: J. Chatterjee and A. N. Das 
Title: A Comparison Of Perturbative Expansions 
IVIanuscript no: BG7510 




FIG. 13 

Authors: J. Chatterjee and A. N. Das 

Title: A Comparison Of Perturbative Expansions 

Manuscript no: BG7510 



